Temperature Dependence of Damping and Frequency Shifts of the Scissors Mode of a 

trapped Bose-Einstein Condensate 
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We have studied the properties of the scissors mode of 
a trapped Bose-Einstein condensate of 87 Rb atoms at finite 
temperature. We measured a significant shift in the frequency 
of the mode below the hydrodynamic limit and a strong de- 
pendence of the damping rate as the temperature increased. 
We compared our damping rate results to recent theoretical 
calculations for other observed collective modes finding a fair 
agreement. From the frequency measurements we deduce the 
moment of inertia of the gas and show that it is quenched 
below the transition point, because of the superfluid nature 
of the condensed gas. 

PACS numbers: 03.75.Fi, 05.30. Jp, 32.80. Pj, 67.90. +z 



The experimental discovery of the scissors mode [0, 
first predicted in a geometrical model has been one 
of the most exciting findings in nuclear physics during the 
last two decades (see 0j for a review). According to the 
geometrical picture, such a mode arises from a counter- 
rotational oscillation of the deformed proton and neutron 
fluids. Extensive studies of this mode in the past two 
decades have investigated the dependence of its strength 
on the nuclear deformation and the relationship with the 
quadrupole collective mode and with the superfluid mo- 
ment of inertia of the nucleus || . 

Recently it has been possible to study the scissors 
mode in Bose-Einstein condensates (BEC) of dilute gases 
H . This transverse mode of excitation has been used to 
demonstrate that the condensate can only flow in an ir- 
rotational manner This work, together with the 

observation of a critical velocity j7j and the occurrence 
of quantized vortices ||p] has provided conclusive evi- 
dence for the superfluidity of a trapped condensate. The 
scissors mode of a trapped BEC has been formally linked 
to the quadrupole operator and to the moment of inertia 
of the superfluid system [jio|| . 

In condensed gases the 'scissors' motion is an irrota- 
tional oscillation of the condensate relative to the static 
trapping potential. The deformation of the system can 
be controlled independently using the trap frequencies 
(unlike nuclei). The most striking difference between the 
scissors mode in nuclei and in condensed gases is that the 
nuclear system is effectively at T = and the scissors ex- 
citation is fragmented because of the interference of the 
orbital with the spin motion, whereas in a trapped BEC 



it is possible to control the temperature of the gas, and 
hence change the fraction No/N of condensed atoms. 

The spectroscopy of the low energy collective modes 
of the condensate has been extended to include measure- 
ments of damping rates and frequency shifts as a func- 
tion of temperature (l^Jl^]. These finite temperature 
measurements test dynamical aspects of the BEC theory 
that require higher order terms in the calculations. As 
yet the results have not been fully explained |^3| . 

In this paper we report a systematic study of the tem- 
perature dependent damping and frequency shifts of the 
scissors mode excitation in a trapped condensed gas. Af- 
ter a brief description of the experimental procedure, we 
present the scissors mode data for the frequency shifts 
and damping rates. We then compare our results with 
the available theoretical calculations for other collective 
modes. Finally we show how the scissors mode frequency 
shifts are related to quenching of the moment of inertia 
of the boson gas. 

We excite the scissors mode using the technique de- 
scribed in our previous paper In summary we pre- 
pare atoms at the desired temperature T in an untilted 
TOP trap {uj x = uj y — 126 Hz, uj z — y/8u x ). We then 
adiabatically tilt the confining potential by a small angle 
4> = 3.4° about the y axis (this reduces ui z by ~ 2%). 
We then suddenly flip the trap angle to —0 to excite the 
scissors mode in the xz plane with an amplitude 8q = 2(f), 
about the new equilibrium position. The condensate os- 
cillates at a single frequency lu sc (Fig. (a)). The low 
temperature limit of the condensate frequency is the hy- 
drodynamic value uj sc — \J <jJ 2 z + u>% ||. The thermal 
component oscillates at two frequencies w+ and u)~, in 
the collisionless limit o;+ = lu z + uj x and oj- = u) z — lo x . 
The high-lying frequency lo + corresponds to an irrota- 
tional quadrupole velocity flow (V x v — 0) that is the 
classical counterpart of the superfluid irrotational oscil- 
lation. The low-lying frequency w_ of the thermal cloud 
is related to rotational flow, and the absence of a simi- 
lar mode for the BEC is the key feature that indicates 
superfluidity. 

To observe the scissors mode oscillation for the thermal 
component we imaged the atomic cloud in the trap using 
destructive absorption imaging. The images of the cloud 
in the trap were fitted with a 2D Gaussian distribution 
from which we extract the cloud angle as a function of 
time (both above and below the critical temperature). 
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To observe the condensate we switch off the trapping 
potential at different times during the scissors mode os- 
cillation, and allow the cloud to expand for 14 ms before 
taking an image. The thermal cloud expands isotropi- 
cally while the BEC component has an elongated (pro- 
late) shape (Fig. |l| (b)). The angle of the expanded 
BEC component is extracted by fitting the two compo- 
nent profile with a 2D double Gaussian distribution jl4j . 
In this way we are able to separate the information on 
the thermal cloud scissors mode and the BEC scissors 
mode. 

The temperature of the ultracold atomic sample is ex- 
tracted from a 2D Gaussian fit to the wings of the ex- 
panded thermal component using the procedure well de- 
scribed in fHf , For temperatures below 0.5To when the 
thermal cloud is not visible, we deduced the temperature 
from the final rf by extrapolating the best fit curve for 
the available temperature vs rf cut frequency data. To 
determine the scaled temperature T/Tq we need to cal- 
culate the critical temperature Tq(N) at each final rf cut 
- in evaporation each rf cut gives a different final num- 
ber of atoms N and hence different To (TV). We calcu- 
lated To (AT) from the measured final number of atoms in 
each time-of-flight picture and the known trap frequen- 
cies, using the formula for a trapped boson gas The 
statistical uncertainty in T/Tq is in the range 3 — 6%, 
but increased to 10% for T/Tq < 0.5 where T was ex- 
trapolated. To check our temperature measurements we 
compared the condensate fraction vs temperature with 
theoretical predictions for both the ideal gas and inter- 
acting gas @ (Fig. (c)). 

The measurements of the frequency and damping of 
the quadrupole modes in previous work [ p~T|JT^ ] relied on 
measurement of the radii of the condensate and the ther- 
mal cloud. These measurements are sensitive to the fit- 
ting model and to shot to shot noise because of their atom 
number dependence. However the scissors mode has very 
well defined initial conditions for the excitation of both 
the BEC and the thermal cloud and the angle oscillations 
of the two components can be measured independently. 
The angles of both condensed and thermal clouds do not 
depend on atom number and could be fitted very reli- 
ably, even close to To to give accurate measurements of 
frequencies and damping. 

Figure H shows three different scissors mode oscillations 
at temperatures of (a) T/T = 1.29, (b) T/T = 0.74 and 
(c) T/Tq — 0.53. The measured frequencies of the ther- 
mal cloud do not change with temperature and agree very 
well with the collisionless frequencies, both above and be- 
low the critical temperature (see also Fig. [|). Both fre- 
quency components have the same amplitude implying 
that energy is shared equally between rotational and ir- 
rotational velocity flow. These thermal cloud excitations 
have very small damping that is independent of temper- 
ature, confirming the collisionless regime of the thermal 
component in our experiment (see also Fig. H). 



The scissors mode oscillation of the BEC component 
occurs at a single frequency u> sc (Fig. || (c) and (b)). 
However, the frequency of oscillation and its damping 
change dramatically with temperature. The low temper- 
ature limit of the frequency was ~ 1% higher than the 
hydrodynamic prediction as shown in Fig. |^. This agrees 
with the frequency we calculated for finite number con- 
densate (2 x 10 4 ) using the method detailed in 

Figure [| shows the experimental data for the temper- 
ature dependent damping rates. Since a finite T theory 
of the scissors mode in a TOP trap is not yet available 
we have compared our experimental results with recent 
theories for other observed collective modes of a BEC. 
Fedichev et al. Jl9| have calculated the Landau damping 
rate for the \m = 2| and low- lying m = modes of BEC 
in the usual TOP trap geometry (u x — u> y , uo z — \^8uj x ) 
for the experimental conditions of the JILA trap |pl| , 
that are close to those of our experiment. In the region 
0.5 < T/Tq < 0.9 their theoretical calculations give a 
damping rate T v for the collective mode v that is pro- 
portional to Tl\x and to the frequency of the excitation 
u) u : 

T v = A^-inoa 3 ) 1 / 2 (1) 
A 4 

where hq is the peak density, fi is the chemical potential 
and A v is a numerical coefficient dependent on the shape 
of the collective excitation. For the quadrupole modes 
they calculated A\ m=2 \ « 7 for the \m = 2| and A m=0 ~ 5 
for the low-lying m=0. The numerical coefficient A v is 
likely to be similar for the scissors mode. Assuming that 
A sc w ^4| m= 2| ~ 7 and scaling the results by Fedichev et 
al. E^] by the frequency factor u> sc /uj\ m= 2\, we obtain the 
theoretical damping rate for the scissors mode (Fig. 
solid line), which is in fair agreement with our data. This 
numerical coefficient for the scissors mode could in prin- 
ciple be evaluated by the method described by Fedichev 
et al. They use an overlap of the wavefunction of 

the mode with the wavefunctions of the thermal exci- 
tations following classical trajectories within the range 
fi < T < 0.9Tb. This method gives a way of calculating 
the damping process for a trapped BEC, the evaluation 
of the trajectories and the stochastic averaging is a nu- 
merical problem of moderate size. Our measurements of 
the damping extend to temperatures where T < M- I n 
this region < T/Tq < 0.5 Fedichev et al. predict damp- 
ing proportional to (T/fi) 3/2 @ but our experimental 
data have a weak dependence on the temperature. 
In a recent paper Rusch et al. (2lJ have used a full sec- 
ond order theory to calculate the damping and frequency 
shift of low energy collective modes for a trapped BEC in 
a spherical geometry . A spherical cloud has no scis- 
sors mode but it is possible to identify this mode with 
the quadrupole mode / = 2,m = ±1. Thus we have 
used their calculations for the I = 2 mode damping rate 
for comparison scaled by the frequency factor lu sc /lui=2 
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(Fig. I], dotted line). These numerical calculations give 
a weakly linear dependence of damping on temperature 
but the experimental data seem to follow a higher power 
law. 

We have not included the observed decrease in damp- 
ing at higher temperatures (above 0.85T ) in the above 
discussion because the uncertainties are large for these 
data points. In addition the frequency shift of the scis- 
sors mode at these temperatures may lead to resonant 
coupling with other modes or with the thermal cloud. 
For example in p0| the lowering of the damping of the 
quadrupole modes close to the transition point, has been 
interpreted as a dynamical coupling between the con- 
densed and thermal component. 

The moment of inertia of a system is defined as linear 
response to a rotational field. Superfluids have a moment 
of inertia less than that of a classical rigid body <d r i g 
with the same density distribution. The scissors mode is 
ideally suited to studying this quenching of the moment 
of inertia below To. Zambelli and Stringari |]l0| linked the 
reduced moment of inertia of the boson gas 6' = 9/ r i g 
to the quadrupole moment Fourier signal Q[oS) for the 
condensate at T = and the thermal cloud at T > Tq. 
Applying their derivation to the scissors mode gives the 
moment of inertia of the condensate and thermal clouds 



as: 



w cond 
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We can use Eq. |2] to deduce the moment of inertia for the 
condensate and thermal cloud separately. The moment of 
inertia for the thermal component is consistent with the 
rigid body value 0J hcrm = 1 for all measured T. For the 
point at 0.99To in Fig. |3| the condensate scissors mode 
has a frequency such that its moment of inertia is close 
to that of a rigid body. Whereas as T — > the conden- 
sate has a quenched irrotational value of the moment of 
inertia characteristic of a superfmid. The lower tempera- 
ture points are in good agreement with the hydrodynamic 
frequency and the reduced moment of inertia assumes 
the superfmid value for our geometry 9^. ond = e 2 = 0.6, 
where e = (w 2 — w 2 )/(o; 2 + ^x) 1S the anisotropy of the 
trap g(J. 

In conclusion we have made precise measurements of 
the frequency and damping of the scissors mode as a 
function of temperature. The observed damping for low 
temperatures is consistent with that expected from mod- 
els of Landau damping, where the damping rate is pro- 
portional to the mode frequency. At temperatures just 
below To the experimental damping rate was less than 
predicted. The frequency shift as a function of temper- 
ature is consistent with a quenching of the moment of 
inertia as the temperature is reduced because of the su- 
perfluidity of the gas. We await a theoretical calculation 
of this shift for comparison. 
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FIG. 2. Angle versus time during the scissors mode, (a) 
Thermal cloud, T/To = 1.29. The scissors mode is character- 
ized by two frequencies of oscillation with equal amplitude, 
(b) Condensate, T/T = 0.74. A heavily damped BEC os- 
cillation is observed at a single frequency, below that of hy- 
drodynamic theory, (c) Condensate, T/Tq = 0.53. Lightly 
damped BEC oscillation at a frequency that agrees well with 
the hydrodynamic value. 
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FIG. 1. (a) At finite temperature, during a scissors mode 
oscillation, the condensate and the thermal cloud oscillate 
at their characteristic frequencies, (b) Typical time of flight 
(TOF) picture during a scissors mode oscillation at finite tem- 
perature: the thermal cloud is isotropically expanded, while 
the condensate has an elongated shape. From a 2D double 
Gaussian distribution fit to the picture we extract the angle 
of the BEC. For this picture N /N = 0.14, T/T = 0.91 and 
6 — —8.5° with respect to the vertical axis, (c) Condensate 
fraction as a function of the reduced temperature T /To. The 
solid line is the interacting theory prediction made by A. Min- 
guzzi et al. for our experimental conditions. The dashed line 
is the ideal trapped gas cubic law prediction. The experi- 
mental points (solid circles) are obtained from TOF images 
as explained in the text. 
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FIG. 4. The temperature dependence of the damping rate 
F normalized to the radial trap frequency uj x for the conden- 
sate (solid triangles) and thermal component (open squares). 
The solid line is the theoretical damping rate of the \m = 2| 
mode calculated by Fedichev et al. (which agrees well with 
the experimental data obtained at JILA for that mode) multi- 
plied by the frequency scaling factor uj S c/^\ m =a\ f° r the scis- 
sors mode. The dotted line is the prediction for the / = 2 
mode by Rusch et al. in a spherical geometry rescaled by the 
frequency factor for the scissors mode. 



FIG. 3. Temperature dependent frequency shift of the scis- 
sors excitation of the condensate (solid circles) and thermal 
component (open circles). Significant negative shifts from 
the hydrodynamic prediction (dashed line) of the condensate 
scissors mode frequency are observed as the temperature in- 
creases. The low temperature frequencies are systematically 
higher than the hydrodynamic value by ~ 1% due to the fi- 
nite number of atoms. The two frequencies of the scissors 
mode of the thermal component do not appear to be temper- 
ature dependent and are in good agreement with collisionless 
predictions (dotted lines). The frequency spectrum has been 
normalized by the radial trap frequency. 
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